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We study the entanglement spectrum of spin- 1/2 XXZ ladders both analytically and numerically. Our analyt- 
ical approach is based on perturbation theory starting either from the limit of strong rung coupling, or from the 
opposite case of dominant coupling along the legs. In the former case we find to leading order that the entangle- 
ment Hamiltonian is also of nearest-neighbor XXZ form although with an in general renormalized anisotropy. 
For the cases of XX and isotropic Heisenberg ladders no such renormalization takes place. In the Heisenberg 
case the second order correction to the entanglement Hamiltonian consists of a renormalization of the near- 
est neighbor coupling plus an unfrustrated next nearest neighbor coupling. In the opposite regime of strong 
coupling along the legs, we point out an interesting connection of the entanglement spectrum to the Lehmann 
representation of single chain spectral functions of operators appearing in the physical Hamiltonian coupling 
the two chains. 

PACS numbers: 75.10.Jm, 03.67.-a 



I. INTRODUCTION 

In the last decade many-body physics has been substan- 
tially enriched by the concept of entanglement, whose exten- 
sive and systematic study originated in the field of quantum 
information theory'. In particular, the notion of the entan- 
glement spectrum^ has led to novel insights in the physics 
of various many-body systems. These include quantum Hall 
monolayers at fractional filling^ quantum Hall bilayers at 
filling factor v — l'^, spin systems of one'^*"'^ and two'^"^' 
spatial dimensions, and topological insulators^^'^^. Other 
topics recently covered encompass rotating Bose-Einstein 
condensates^"*, coupled Tomonaga-Luttinger liquids""", and 
systems of Bose-Hubbard"'' and complex paired superfluids"^. 

In Ref. 17 Poilblanc observed that chain-chain entangle- 
ment spectra in two-leg spin ladders are remarkably similar 
to the energy spectrum of a single Heisenberg chain. Fur- 
thermore he found that the fitted effective inverse temperature 
depends on the ratio of the leg to the rung couplings and van- 
ishes in the limit of strong rung coupling. In a parallel study 
of entanglement spectra of quantum hall bilayers at = 1 one 
of the present authors observed a similarly striking analogy 
between the entanglement spectrum of a single layer and the 
energy spectrum of a single physical layer siv ~ \/2}^ 

In this paper we study the entanglement spectrum of two 
coupled XXZ chains analytically in two limiting cases a) the 
case of strong rung coupling and b) the case of weak rung 
coupling. In case a) we find that the entanglement spectrum 
is described to leading order by an entanglement Hamilto- 
nian of the XXZ form, albeit with an in general renormal- 
ized effective anisotropy (and thus independently reproduc- 
ing the recently posted result of Ref. 18). For the particular 
cases of XX and Heisenberg ladders the anisotropy is unal- 
tered, and we arrive at the conclusion that in these cases the 
entanglement Hamiltonian is indeed the physical Hamiltonian 
restricted to a single chain. Moreover, we derive explicit re- 
sults for the next-to-leading order exhibiting deviations of the 
entanglement Hamiltonian from nearest-neighbors XXZ cou- 
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FIG. 1 : (Color online) Two leg spin ladder consider in this paper The 
chains are labeled hy v = 1,2, while the runs are labeled by m = 
1, . . . L. We consider the entanglement spectrum in the illustrated 
chain-chain bipartition. 



pling. In the second case (b) we point out an interesting con- 
nection between the entanglement spectrum and the Lehmann 
representation of single chain spectral functions of the oper- 
ators contained in the physical Hamiltonian coupling the two 
chains. 

We consider an XXZ spin ladder Hamiltonian T-L = Ho + 
Hi with 



Ho — «/i 



rung 



,(1) 



describing the coupling along the rungs and legs, respectively. 
The sites are labeled by (to, v), where to S {1, . . . , L} de- 
notes the position within chain = 1,2. All spin-1/2 op- 
erators are taken to be dimensionless such that Jrung^ -^log 
have dimension of energy, and A is a dimensionless XXZ 
anisotropy parameter. In the following we will always assume 
a ladder of length L with periodic boundary conditions if not 
specified otherwise. In the following we limit our discussion 
to the antiferromagnetic regime Jiog, Jrung > 0. Apart from 
the special case Jicg = 1, Jrung — (corresponding to two 
decoupled critical S = 1/2 chains), the system is gapped in 
this regime, with exponentially decaying correlations. 
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to their quantum number^^. In the limit J\cg^/ J^^ 



0+ 



all the 2^ levels of the entanglement spectrum collapse onto 
a single value Lin 2, and for finite J\cg/ J-mng they start to 
spread. For Jiog/ Jrung It 1 the spectrum reaiTanges, and a no- 
table feature is that the lowest levels above the singlet ground 
states are a sizable set of = 1 triplets, displaying a common 
slope of 2 when plotted as a function of ln( Ji^g/ Jrung)- Note 
that this is in contrast with the behavior at small Jiog/ Jrung, 
where there is only one triplet located below the second sin- 
glet in the entanglement spectrum. 

In the lower panels of Fig. 2 we highlight the behavior in 
the two limit cases by appropriately rescaling and shifting the 
entanglement spectra. In the lower left panel we also show 
the energy spectrum of a L = 8 Heisenberg chain (multiplied 
by two) to demonstrate the remarkable agreement between the 
entanglement spectrum and the energy spectrum in this partic- 
ular limit. As just discussed the structure of the entanglement 
spectrum is somewhat different in the opposite limit, and the 
analogy with the energy spectrum of a single chain apparently 
breaks down. 

We now provide an analytical justification of the reported 
behavior based on perturbation theory around the two limits 
(a) and (b) described above. 



III. STRONG RUNG COUPLING LIMIT 

Let us now treat Hi as a perturbation to "Ho with antiferro- 
magnetic coupling, Jrung > 0. The unperturbed ground state 
reads 



FIG. 2: (Color online) Overview of the entanglement spectrum of 
the reduced density matrix of a single chain (cf. setup in Fig. 1) in 
a Heisenberg (A = 1) spin ladder for different Si^t = 'S'tot sec- 
tors as a function of Jieg/Jrung- The system size is L = 8. While 
in the upper panel the bare entanglement spectrum is shown, in the 
lower panels the leading asymptotic behavior in the two limiting 
cases Jicg/ Jrung — >■ and Jrung /Jicg — > is highlighted by appro- 
priate shifts and rescaling indicated on the y— axes. The filled dots in 
all panels denote numerical entanglement spectrum levels, while the 
constant lines in the lower left panel highlight twice the eigenvalues 
of a single L = 8 S = 1/2 Heisenberg chain. The grey region in 
the lower right panel indicates the part of the rescaled entanglement 
spectrum affected by the finite (double) precision arithmetic used in 
the numerical calculations. 
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using obvious notation for singlet and triplet states on the 
rungs. 
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II. EXACT DIAGONALIZATION RESULTS 



A. First order 



Before discussing our analytical results in the two limiting 
cases we present an overview of the full Jicg /Jrung depen- 
dence of the entanglement spectrum ^ — log(A) obtained 
from the eigenvalues A of the reduced density matrix of a sin- 
gle chain in a Heisenberg (A = 1) spin ladder (cf bipartition 
setup shown in Fig. 1). In the region Jicg/ Jrung < 1 selected 
numerical spectra have aheady been presented in Refs. 17,20. 

In the upper panel of Fig. 2 the bare entanglement spectrum 
for a L = 8 ladder is shown as a function of Jiog/ Jrung (note 
the logarithmic x-axis), where the levels are labeled according 



The first-order correction to the ground state can be ob- 
tained by elementary calculation. 
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where the dots denote singlet states on each rang not expUcitly 
specified. The reduced density operator is obtained by tracing 
out one of the legs from p — (|0) + 11))((0| + (1|) and is given 
within first order in Jicg/ Jrung by 




|^('5'm'5'm+l + h.C.) 



m+1 



(8) 



with L being the number of rungs. Again within first order 
perturbation theory, this result can be formulated as 



1 



/'red = ^ exp 



with Z = Trexp(— 'H^J,!.) being a partition function. The 
entanglement Hamiltonian "Hpni is given as: 
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and is of nearest neighbor XXZ form, with a renormalized 
anisotropy parameter. 



A = i(A + A2) 



(11) 



Note that A = (XX case) and A = 1 (Heisenberg interac- 
tions) are invariant, i.e. the entanglement Hamiltonian is pro- 
portional to the physical Hamiltonian restricted to the block, 
as observed numerically for A = 1 earlier in Ref. 17. Since in 
these two specific cases the physical Hamiltonian on a chain 
and the entanglement Hamiltonian are simply proportional to 
each other (to first order), one can define an ad hoc inverse 
temperature 
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which can, within the same order, be reformulated as 
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with Z — trcxp(— "Hcni) and 
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Thus, also within second order Jiog/ Jrung the entanglement 
Hamiltonian is quite similar to the physical Hamiltonian on 
a single chain. The second order correction contains a renor- 
malization of the amplitude of the nearest neighbor coupling 
and the appearance of an unfrustrated ferromagnetic next 
nearest neighbor Heisenberg coupling. The latter observation 
is in accordance with the numerical findings in Refs. 18,20. 
We expect this picture to hold also at higher order According 
to general linked cluster ideas, sites at distance r will only 
be able to couple starting at order r with a leading ampli- 
tude proportional to ( -p£fi- j . Given the Heisenberg nature 



of the leg interactions, to leading order a Heisenberg interac- 
tion J2m ^mSm+r IS the Only symmetry-allowed interaction. 

In Fig. 3 we present a comparison between the rescaled 
numerical entanglement spectra for L — ^ and the perturba- 
tive result in the form of the energy spectrum of Eq. (15) for 
Jicg/ Jrung < 1- The second order corrections only seem to 
become visible for this system size for Jieg/ Jrung 0.1, and 
the initial deviations from a straight line in the numerics are 
correctly described by the second order corrections. However 
for larger Jicg /Jrung ratios further neighbor Heisenberg cou- 
plings and multi-spin interactions become non-negligible, as 
observed in Ref. 20. 



IV. LIMIT OF WEAKLY COUPLED CHAINS 



B. Second order 



We now discuss the opposite limit of weakly coupled 
chains. At the starting point Jmng/ Jeg = the ground state 
is the product of the ground states of the individual chains. 



The second-order contribution to the ground state is some- 
what lengthy and given in appendix A. For general anisotropy 
the result does not seem to be amenable to a simple interpre- 
tation. In the isotropic case A = 1, however, one finds up to 



IV'o) = |0)a|0)j 



(16) 



In this limit the entanglement spectrum is trivial, composed 
of only a single value ^ = 0, 
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FIG. 3: (Color online) Filled circles: numerical entanglement spec- 
tra obtained for a L = 8 system. For each Jicg/Jrung value the 
entanglement spectrum has been rescaled to lie in the interval [0, 1]. 
Empty circles: rescaled energy spectrum of the analytical second or- 
der entanglement Hamiltonian Eq. (15). 



To first order the wave function reads: 
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FIG. 4: (Color online) Filled diamonds: entanglement spectrum 
prediction for a L = 8 Heisenberg spin ladder, based on the first 
order perturbation theory in the coupling between the two chains 
[Eqn. (22)]. Empty diamonds: S = and S = 1 entanglement 
levels obtained from exact diagonalization at Jrung/ Jicg ~ 0.0001, 
cf. upper panel of Fig. 2. The grey region denotes the region where 
the numerical entanglement spectrum is affected by the finite preci- 
sion arithmetics and where no numerical data is shown. 



|V'i> = [l-i?Qo^o]|^/'o), 



(17) 



where Hq is the Hamiltonian coupling the two chains [cf. 
Eq. (1)], Qq is the projector onto the subspace orthogonal 
to \ipo), and R denotes the resolvent operator R = {Hi — 

After a straightforward calculation one obtains the follow- 
ing expression for the first order wavefunction lipi) in the 
product basis of chain eigenfunctions: 
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where k runs over the lattice momenta of a single chain, a runs 
over the spin components {x,y}, n and n' label the eigen- 
functions of the two isolated chains, while A, 



in') 



denotes the 



single chain excitation energies — Eq. One recognizes 

that the nontrivial part of the wave function is composed of 
contributions which also enter the Lehmann representation of 
the 5""(fc, Lj) spectral functions of a single chain, 

S''{k,uj) = ^|c^-|2x,5(l.-A„), (20) 

n 

ct-" = {n\Sm, (21) 



where the matrix elements cji'" enter as 
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Given the wave function in this form, the entanglement 
spectrum can simply be obtained by a singular value decom- 
position of the matrix ]«,«', which amounts to finding the 
Schmidt decomposition in this bipartition. The singular val- 
ues are thus the Schmidt values, and when squared they cor- 
respond to the eigenvalues of the reduced density matrix in 
the same bipartion setup. Due to the translational invariance 
along the chains as well as the total preserving form of 
the XXZ Hamiltonians, the wave function matrix exhibits a 
block diagonal structure in k and only n and n' sectors with 
total = 0, ±1 appear. In the specific case of Heisenberg 
Hamiltonians (A = 1) one can then infer by virtue of the 
Wigner-Eckart theorem that only total spin S — 1 entangle- 
ment levels appear above the singlet ground state to leading 
order. In addition the singlet ground state leads to vanishing 
spectral weight at fc = and therefore the absence of 5 = 1 
entanglement levels at fc = to the expansion order consid- 
ered. 

In Fig. 4 we show a comparison between the numerical en- 
tanglement spectrum for i = 8 at Jmng/ J\cg = 0.0001 and 
the prediction based on Eqn. (22). The filled symbols rep- 
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resent all the predicted entanglement levels appearing at this 
order. The numerical results match perfectly the analytical 
prediction up to the finite precision threshold of the numerics. 

V. DISCUSSION AND CONCLUSION 

In summary, we have investigated the entanglement spectra 
of spin- 1/2 XXZ ladders using both exact numerical diago- 
nalizations and perturbation theory approaches. In the limit 
of strong rung coupling perturbation theory predicts in lead- 
ing order an entanglement Hamiltonian being also of nearest- 
neighbor XXZ form with a renormalized anisotropy parame- 
ter. For XX and isotropic Heisenberg ladders no such renor- 
malization takes place and in this case one can define an ad 
hoc effective temperature proportional to the Jrung/ J\ce ratio. 
Moreover, in next-to-leading order the entanglement Hamilto- 
nian exhibits spin couplings of longer range, and all the above 
findings are in perfect agreement with our numerical exact di- 
agonalization results. The first order perturbation theory re- 
sult has also been obtained in parallel work by Peschel and 
Chung"*. Similar numerical observations have already made 
for quantum Hall bilayers at filling factor u — l'-. It remains 
an open task to devise analogous perturbational arguments for 
such systems with long-ranged interactions . 



In the opposite regime of strong coupling along the legs, we 
point out an interesting connection of the entanglement spec- 
trum to the Lehmann representation of single chain spectral 
functions of the operators contained in the physical Hamilto- 
nian coupling the two chains. This aspect also holds true for 
other ladder systems, such as e.g. Hubbard ladders, where the 
entanglement spectrum will be determined by elements of the 
single particle addition and removal spectral functions in the 
weakly coupled chain regime. 
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Appendix A: Second-order correction to tlie ground state and 
the reduced density matrix 

The second-order correction to the ground state can be ob- 
tained as 
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Note that no terms with an odd number of triplets occur In the expressions, 
isotropic Heisenberg case A = 1, the second-order contribu- 
tion to the reduced density matrix is the sum of the following trneg (|0)(2| + |2)(0|) 




Now using the identities for spin-l/2-operators, one derives the result (13). 

^-S.Ss = (5^16^2) (525^3) + (5^25^3) {S1S2) (A4) 
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and 

3 1 



- ^SiS2 = (^3182) (A5) 
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